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We investigate Ua (1) breaking above T c in terms of the Dirac spectrum on configurations with 
(2+l)-flavors, using the HISQ action. The strange quark mass is at its physical value. We use 
several light quark masses corresponding to the Goldstone pion masses in the range of about 
115 - 230 MeV on lattices of size 32 3 x8 and 48 3 x8. We calculate the 100 lowest-lying Dirac 
eigenvalues at temperatures below and above T c . We investigate the volume dependence of the 
Dirac eigenvalue density to determine whether there is a gap around zero, which can appear if 
Ua{1) symmetry is restored in the chiral symmetric phase. We also investigate the quark mass 
dependence of the Dirac eigenvalue density at zero and check whether there is a linear behavior 
that would signal the Ua(1) breaking above T c . 
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1. Introduction 

Quantum chromodynamics (QCD) with Nf flavors is invariant under a global flavor symmetry 
of SUi(Nf) x SUit(Nf) x t/y(l) x Ua(1) in the limit of vanishing quark mass. It is also known 
that the SUi(Nf) x St/«(iV/) chiral symmetry is spontaneously broken in the vacuum and the Ua(1) 
symmetry is explicitly broken due to the axial anomaly. However the chiral symmetry is expected 
to be restored at high temperature. On the other hand, although the anomaly effect is present at 
any finite temperature, it is suppressed at high temperature and accordingly the Ua(1) symmetry is 
effectively restored. Here an important question is when and how such an effective Ua (1) symmetry 
restoration occurs since it can influence the order of the chiral phase transition. In the case of Nf = 2 
it has been suggested ^ that if the C/a(1) symmetry remains broken at the critical temperature T c 
of the chiral phase transition, the phase transition is of second order belonging to a universality 
class of a 3-dimensional 0(4) spin model. On the other hand, if both the chiral and the Ua(1) 
symmetries are restored at the same time, the phase transition should be of first order. 

To study the QCD phase structure and properties of strongly interacting matter, where non- 
pertabative effects play an important role, lattice QCD simulations are a suitable way. Some lattice 
QCD studies with (2+l)-flavor staggered type quarks |3j] have indicated at least 0(N) scaling in 
the analysis of the magnetic equation of state, although only the 0(2) symmetry is preserved. Sim- 
ilarly in the case of 2-flavor 0(a) improved Wilson quarks [Q, 0(4) scaling has been reported. 
This suggests the second order phase transition scenario in the case of Nf = 2 massless quarks. In 
our previous report ^ we investigated the temperature dependence of the Dirac eigenvalue distri- 
bution and discussed a scaling behavior in the (2+l)-flavor case with the highly improved staggered 
quark (HISQ) action There we found that the eigenvalue density around zero, which is related 
to the chiral order parameter, has a crossover behavior and a power law dependence similar to 
what is expected from the critical exponent of the 0(2) or 0(4) universality classes. Moreover, we 
found no clear evidence for a gap around zero in the Dirac eigenvalue density at high temperature, 
which supports Ua(1) breaking above T c . Recently the restoration of the Ua(1) symmetry in terms 
of the Dirac eigenvalue distribution has been discussed also with domain wall and overlap [§] 
fermions, which have exact chiral symmetry on the lattice. With domain wall fermions there re- 
main small eigenvalues even at high temperature. Whereas with overlap fermions a gap around 
zero seems to appear above a certain temperature. However, the lattice sizes in these studies are 
not large enough for definite conclusions and more detailed investigations are needed. 

In this report we continue our study on the Dirac eigenvalue distribution at finite temperature 
with the same setup as in ^ and give more details about its volume and quark mass dependences 
and the effective £7^(1) restoration. In the next section we discuss details of our study. Then we 
show some numerical results in Sec. ^] and finally give conclusions in Sec. ||. 

2. Effective Ua(1) restoration and the Dirac eigenvalue density 

To discuss the effective restoration of the Ua(1) symmetry, we consider the chiral condensate 
(yy) and several susceptibilities % defined in the 4-dimensional Euclidean space-time as follows: 

T dlnZ T , i. 
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Xs = Xcon = —(TviM-'M- 1 )), (2.2) 

Zd lsc = \ [((TrM- 1 ) 2 ) - (THM- 1 )} 2 ] , (2.3) 
_ T d 2 lnZ _ d{\jf\jf) _ 

Xo~ — T7 -i 9 — -i — Xcon T" ^disc j (2-4) 

V dmr am 

^ = ^(Tr(M-SM-V 5 )), (2-5) 

where Z is the QCD partition function and V, T, m and M are volume, temperature, quark mass 
and the fermion matrix, respectively. If the chiral symmetry is unbroken, (yy) = and Xn = Xa, 
and accordingly %%—%s = Zdisc- In addition, if the E/a(1) symmetry is also unbroken, % n = %s an d 
Xk~Xs = Zdisc = 0. Here {i/Zi//) and % % — %§ can be given in terms of the Dirac eigenvalue density 
p(A,m) = (T/^^S^^-A)) as follows: 

r°°2mp(X,m) , 
Jo A z + m z 

r°4ra 2 p(A,m) . 

X*-Xs= jrrz 2^2 dX ' V n) 
Jo (A z +m z ) z 

where we take V — > °°. Then, if m — > is also taken, the Banks-Casher relation [[K]] (v/'i/f) = 



7ip(0,0) follows from (2J5), which means p(0,0) = in the chiral symmetric phase. Thus eigen- 
values close to A = should be suppressed around T c even if the quark mass is finite. On the other 
hand, ( p7| ) tells us that whether % n — %s vanishes in the chiral limit or not depends not only on 
p(0,0) = but also on the approach of p(A,m) to the origin in the limits of m — > and A — > 0. 



Since only eigenvalues around A = can contribute to the integral (2.7) when m —> is taken, one 
possibility to let % n — Xs become zero is that p (A , m) has a gap around A = 0. 

Here, to consider other possibilities, we assume that the dominant contribution to p (A , m) at 
small eigenvalues and quark masses is of the form cm a -Ptf. In this case the above two quantities 
are approximately rewritten as 

{yy} ~2cm a h(p), %n~ Xs ^cm"" 1 ^), (2.8) 

with 

r A/m X P 
/() 

In this study we examine whether the eigenvalue distribution p(A,m) remains finite at A = 0, i.e. 
we consider the case (5=0. The dominant low eigenvalue structure of p(0,m) thus is assumed to 



W = L TTTZaskdx. (2.9) 



ISC 



be p(0,m) = cm a . It is evident from (2JS) that we need a > 1 in order to have % K — %§ = Xd 
finite above the critical temperature. Moreover, only for a = 1 does Xn~Xs = Zdisc stay non-zero 
and thus would provide a sufficient signature for ( 1) symmetry breaking. Of course, a > 1 or 
c = are viable possibilities. In that case any remnant £/a(1) symmetry breaking would have to 
arise either from the structure of p (A , m) at A > or a possible zero mode contribution p (A , m) ~ 
m 2 8(X). 
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Table 1: The parameters and statistics of the numerical simulations. j3, a, T, m s , m\, N s and N T are the lattice 
gauge coupling, lattice spacing, temperature, bare strange quark mass, bare light quark mass, spacial lattice 
size and temporal lattice size, respectively. 





a [fm] 


T [MeV] 


m s a # confi 


gurations 










mi/m s = 1/10 


1/20 


1/27 


1/40 








x Ni = 32 3 x 8 








6.195 


0.1932 


127.7 


0.0880 


702 


_ 


360 


6.245 


0.1832 


134.6 


0.0830 


792 




822 


6.260 


0.1803 


136.8 


0.0810 


768 


- 


834 


6.285 


0.1756 


140.4 


0.0790 


714 




882 


6.315 


0.1702 


144.9 


0.0760 


684 




936 


6.341 


0.1656 


148.9 


0.0740 


612 




402 


6.354 


0.1634 


151.0 


0.0728 


534 


1158 


600 


6.390 


0.1573 


156.7 


0.0694 


600 


1458 


600 


6.423 


0.1520 


162.2 


0.0670 312 


696 


1362 


540 


6.445 


0.1486 


165.9 


0.0652 


600 


1146 


522 


6.460 


0.1463 


168.5 


0.0640 


582 






6.488 


0.1422 


173.4 


0.0620 


678 






6.550 


0.1335 


184.8 


0.0582 


198 






7.150 


0.0747 


330.1 


0.0320 


594 






48 3 x8 


6.423 


0.1520 


162.2 


0.0670 


144 






6.445 


0.1486 


165.9 


0.0652 


156 






7.150 


0.0747 


330.1 


0.0320 


306 







3. Numerical results 

We use the tree level improved gauge action and the HISQ action [|7|] which can reduce cutoff 
effects due to the taste-symmetry violation better than any other staggered type quark action. Part of 



our gauge configurations have been generated by the HotQCD collaboration [11]. We perform our 
simulations mainly on 32 x 8 lattices and for a few temperatures on 48 3 x 8 lattices to investigate 
the volume dependence. The scale was set from the kaon decay constant fx- Our strange quark 
mass m s is fixed to its physical value and we have several light quark masses mi which have values 
of m s /40, ra.v/27, m s /20 and m s /10 corresponding to the lightest (Goldstone) pion of about 115, 
140, 160 and 230 MeV, respectively. Here to/ = m s /21 gives an almost physical pion mass and 
the corresponding pseudocritical temperature in the continuum limit has been estimated at 154(9) 
MeV by using 0(4) scaling fits to the chiral condensate and susceptibility. The details of the 
determination of the lattice spacing, the strange quark mass and the pseudocritical temperature 



have been discussed in Ref. [11]. Each 10th trajectory is chosen for measurements after skipping 



at least 500 trajectories for thermalization. Statistical errors are estimated by the jackknife method. 
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Our simulation parameters and statistics are summarized in Table |TJ 

We calculated the lowest 100 positive eigenvalues of the staggered (HISQ) Dirac operator for 
all parameter sets since the staggered Dirac operator is anti-hermitian and has only pure imaginary 
eigenvalues in complex conjugate pairs. We evaluate the eigenvalue density by binning eigenvalues 
in some small intervals for each configuration and normalize it so that / dXp{X,m) = (T/V)x# 
eigenvalues. 

In our previous report [|]] we found suppression of p ( A , m) around A = and lack of eigen- 
values very close to zero above T c . At the same time, there was also a tail approaching the origin 
in p(A,m) even at our highest temperature, i.e. 330.1 MeV. Since small eigenvalues could suffer 
from finite volume as well as topological effects which should vanish in the thermodynamic limit, 
investigating the volume dependence of p(A,m) is needed to judge whether there is a real gap 
around A = in p(A,m). In the left side of Fig. |l] the volume dependence of p(A,wi) at T = 330.1 
MeV is shown; p(A,m) seems to have a quite small volume dependence in a region of A < 0.05, 
although statistics is not enough to see it clearly. We also check for the effect from the (would-be) 
zero modes of topological non-trivial configurations whose contribution to the eigenvalue density 
at zero would vanish in the thermodynamic limit. The right side of Fig. [l] shows the same result 
as on the left side but for topological trivial configurations. Here the topological charge Q is de- 
termined form a lattice FF operator with HYP fat links It is clearly found that contribution 
form small eigenvalues is more strongly suppressed and eigenvalues less than A ~ 0.03 disappears 
in this case. However, since p (A , m) seems to decrease continuously and exponentially, it is possi- 
ble that the tail of p(A,m) continues below the lowest value of our p(A,m), when a more precise 
determination is made. Increased statistics and bigger volumes are needed to make more definite 
statements. 

As discussed in the previous section, the quark mass dependence of p (0, m) is also important to 
understand C/a(1) restoration above T t . To estimate p(0,m) we fit p(A,m) to a cubic polynomial. 
In this fit analysis we choose the same bin size of 0.0005 for all eigenvalue densities which we 
calculated in order to control a systematic uncertainty coming from the bin size dependence. The 
lower bound of the fit range is fixed to the lowest bin and the upper bound is chosen so that # 2 /dof 
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Figure 2: The temperature dependence (left) and the quark mass dependence (right) of the Dirac eigenvalue 
density at A = 0. 



is as close to unity as possible. Then our # 2 /dof is in the range from 0.88 to 2.1. In the left side of 
Fig. ||the temperature dependence of p(0,ra) is shown. It can be seen that p(0,m) monotonically 
decreases as the temperature increases and, within statistical errors, it seems to vanish around 
T = 160 - 170 MeV at least for mi/m s = 1/40, 1/27 and 1/20. Here we note that even above 
the crossover transition temperature T c , as determined from the chiral condensate or susceptibility, 
the condensate will stay non-zero for non-zero quark masses. The eigenvalue density p(0,m) thus 
may stay non-zero and approaches zero smoothly also above T c . The right side of Fig. || shows 
the quark mass dependence of p(0,m) at four temperatures around T c , namely T = 151.0, 156.7, 
162.2 and 165.9 MeV. At first, p(0,m) seems to approach zero in the chiral limit at the lower 
two temperatures as expected in the chiral symmetric phase. However it is difficult to determine 
whether p (0, m) approaches zero linearly so far. To do so we need simulations with smaller quark 
mass. On the other hand, data for the higher two temperatures are very close to zero already at 
mi/m s = 1/40. However, data for two different volumes at mi/m s = 1/20 show finite volume 
effects which indicate an increase of the density in the thermodynamic limit. This means that our 
results at T = 166.8 and 170.5 MeV underestimate the infinite volume result. We thus need larger 
lattices to get more accurate results. In addition, the eigenvalue density around zero at temperatures 
close to T c is sensitive to fit Ansatze since it changes quickly with respect to X. Hence using any 
other fit Ansatze will help to estimate systematic uncertainty. 

4. Conclusions 

We study Ua ( 1 ) symmetry breaking in the chiral symmetric phase in terms of the Dirac eigen- 
value distribution. In our (2+l)-flavor dynamical simulations, using the tree level improved gauge 
action and the HISQ action, we investigate the volume and quark mass dependence of the eigen- 
value density. 

Even in the eigenvalue density at T = 330.1 MeV, we find that there is a tail approaching the 
origin with only a small volume dependence. Moreover the lower end of the tail vanishes when 
only topological trivial configurations are used. However it is not clear whether there is a gap or 
whether a tail exists and would be visible with more statistics. We also find that the eigenvalue 
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density at zero seems to go to zero in the chiral limit as expected for an unbroken chiral symmetry, 
at T = 151.0 and 156.7 MeV, which are temperatures around T c . We also try to see, if the eigenvalue 
density at zero has a linear quark-mass dependence, which would be sufficient to signal breaking 
of the C/a(1) symmetry in the chiral symmetric phase. However it is difficult to identify it so far. 

Increasing statistics and simulations with larger volumes and smaller quark masses are need 
to make the above points clear. 
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